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Abstract 

We generalize the Brundobler-Elser hypothesis in the multistate Landau-Zener problem to the 
case when instead of a state with the highest slope of the diabatic energy level there is a band 
of states with an arbitrary number of parallel levels having the same slope. We argue that the 
probabilities of counterintuitive transitions among such states are exactly zero. 
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The multistate Landau-Zener problem has been an active field of research during the last 
decade with various applications in condensed matter and atomic physics. The two state 
problem with linear time- dependence of diagonal elements of the Hamiltonian was solved 
exactly by Landau and Zeneri^. The simplest generalization of the two state problem is the 
Schrodinger equation of the form 

iip(t) = (A + Bt)i>(t) (1) 

where A and B are Hermitian N x N matrices with constant elements. The Hamiltonian of 
the model is H = A + Bt. The matrix B can be always chosen diagonal. The goal of the 
theory is to find the transition probabilities, namely the squared elements of the scattering 
matrix lim \Sij(t',t)\ 2 where % and j enumerate eigenstates (the so called diabatic 

f— »+oo ,t—*— oo 

states) of the matrix B with time-dependent diabatic energies E t (t) = But + An = (3{t + Oj. 
Nondiagonal elements of the matrix A that couple diabatic states with the same slopes $ 
can always be made zero by a change of the basis. It is convinient to visualize the time- 
dependence of diagonal elements of any such model in the time-energy diagram like the one in 
Fig^ When all crossing points are well separated one can try to solve the problem naively by 
a successive application of the two state Landau-Zener formula at every two level intersection. 
Even in this approximation the dependence of transition probabilities on parameters can be 
very complicated since amplitudes of different paths leading to the same final states can 
interfere. The task becomes even more complicated when more than two levels can be close 
to each other simultaneously. Then even approximate estimates become very sophisticated^. 
In spite of this complexity, there have been a number of remarkable efforts to solve the model 

exactly at least for some special choices of parameters. Generally this requires nontrivial 
approaches because to solve the n-state model one must consider a n-th order differential 
equation with time-dependent coefficients. 

Although a few important classes of exactly solvable models of the type (JTJ have been 
known for a long time^ the interest toward exact results in the multistate Landau-Zener 
problem has grown up after the work of Brundobler and Elser 6 , who noticed that for any 
model of the form (jlj there are elements of the transition probability matrix that can be 
found by a simple application of the two state Landau-Zener formula at every intersection of 
diabatic energies. Particularly, they presented an empirical formula for the diagonal element 
of the scattering matrix for the state whose diabatic energy level has the highest slope, i.e. 
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FIG. 1: Diabatic energies of a 5-state Landau-Zener model. The choice of parameters is as follows, 
P 1 = j3 2 = f3 3 = l, 4 = 0, fa = -0.8, a\ = 0, a 2 = 0.3, a 3 = 0.5, a 4 = 0, a 5 = 0.4. 

if k is the index of the state with (3k = max(/3i . . . (3n) or (3k = min(/?i . . . (3n) then 



The formula (J2J) is confirmed by all known exactly solvable models with finite number of 
state a 4 i 5 ' 7 i 8 i 9 i 10 and by multiple numerical checks. The authors of& speculated that this find- 
ing probably indicates that the whole problem can be solved exactly or at least can be 
understood in terms of the two-level crossings. Various exact solutions and approximations 
seem support this ide a lc>l11 . Recent work — demonstrated that (J2J follows from a simple 
analytical continuation of the asymptotic solution into the complex time, though such a 
procedure fails to predict correctly other elements of the scattering matrix. The goal of the 
present work is to demonstrate that the Brundobler-Elser hypothesis can be generalized to 
some nondiagonal elements of the scattering matrix and to explain why analytical continu- 
ation of amplitudes into the complex times provides correct predictions for some elements 
of the scattering matrix. 

Assume that instead of one state with the highest slope of diabatic energy level there is a 
band of an arbitrary number of states having the same highest slope so that diabatic energies 
in this band are different only by constant parameters a m . If we assume a " semiclassical" 
approximation where transition between any two states happen only at the corresponding 
crossing point of their diabatic energies then there are elements of the transition probability 
matrix that would be zero in this approximation. Such transitions, if happen, are called 
counterintuitive transitions 13 . Thus, in the model shown in Fig.l, transitions from the state 
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1 to states 2 and 3 and from the state 2 to the state 3 are counterintuitive. 

Generally for the model (JTJ) , if p m = P n = max(/3i . . . Pn) then the transition from the 
state m to the state of the same band n would be counterintuitive if a m < a n . Correspond- 
ingly, if (3 m = f3 n = min(/3i . . . Pn) then the transition is counterintuitive if a m > a n . We 
argue that in the multistate Landau-Zener model with linear time-dependence of diabatic 
energies such counterintuitive transitions have exactly zero probability i.e. without assum- 
ing any semiclassical approximation for any model of the type (JIJ), if the transition from the 
state m to the state n is counterintuitive, then 

\S nm (+oo, -oo)| = (3) 

The "no-go" formula Q and the Brundobler-Elser conjecture (J2J) can be understood by 
the approach similar to the one used by Landau in the two state calculations 1 . Since we are 
interested in the asymptotic magnitude of the amplitudes we can analytically extend the 
evolution (JJJ to imaginary time and choose the evolution path so that always \t\ — > oo. The 
distances between instantaneous eigenenergies 6j(t) of the Hamiltonian remain always large 
in this case, namely of the order of — Pj)t\ » \Aij\ for the states i ^ j and hence we 
can use the adiabatic approximation 

A(t) = e- iS * t o eim ^(t ) (4) 

where the state ipi has the leading asymptotic ipi ~ exp(— i(3it 2 /2) at t — > — oo. 

The approximation (JIJ becomes exact in the limit t — > oo but it is valid generally only if 
there are no other solutions that become exponentially large in comparison with the state ipi 
to which it is applied. Suppose that the state ipo has the largest slope of the diabatic energy 
Po at t —>■ — oo and is initially occupied. In this case it is convenient to choose the time-path 
as shown in Fig|2]with t = Rexp(i(p) where R — > oo and <fi decreases from ir to zero. One 
can always change variables so that /3 = and $ < for states with slopes /3j ^ (3<jL. When 
(p changes in the interval from 3tc/A to 7r/4, the amplitudes of states with slopes f3i < 
are decreasing exponentially and become suppressed by the factor exp(C(0)/3j|t| 2 /2) where 
Pi < and C(<p) is a positive coefficient that depends only on the angle. We choose the 
asymptotics so that at the angle <fi = 3ir/4 the state i/jq is dominating over all others, i.e. is 
exponentially large in comparison to them. Then the states with $ < should not affect 
the adiabatic approximation in the interval 37r/4 > <fi > it/4 since they can only decrease 
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FIG. 2: The deformed time contour for the evolution from large negative to large positive times 
with t = Rexp(icft), R — ► oo, < <j) < 7r. 

there. One can see that the condition that at <fi = 37r/4 the state ipo is dominating also leads 
to the vanishing of the amplitudes of other states with fa < in the interval it < <p < 3n/A 
so that it is not forbidden to choose \ipo(— oo) | = 1 and \i>i{— oo)| — > 0, (i ^ 0). 

At the last part of the contour 7r/4 > > amplitudes of states with fa < grow from 
almost zero value, but at <p = time becomes real and hence amplitudes cannot be larger 
than unity. So in this part of the contour such amplitudes have not enough time to become 
exponentially large. It means that they still remain small or comparable with ^ at this 
interval and the formula (J3J) should be valid for the state ipo during the whole evolution. 
Substituting the energy up to the first order correction in l/\t\ 

IA:nl 2 



e (t) ~ a + 



(fa ~ fa)t 



(5) 



into the formula for the transition probability 

\M+oo)\ 2 



I'S'ool 2 



exp —2Im( / e (t)dt) 



(6) 
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we find the Brundobler-Elser result (0). It is clear from this analysis why the formula (0) is 
generally not valid for other diagonal elements of the scattering matrix. If an initially filled 
state does not have the highest slope of the energy level there are states with higher slopes 
that grow exponentially and become large in the interval 37r/4>0>7r/4of the contour so 
that the adiabatic approximation becomes invalid in application to ipQ. To treat this case 
properly, one should investigate the Stokes phenomenon near all crossing points of diabatic 
energies^. 



This analysis becomes more complicated if there is more than one state having the same 
largest energy slope (3q. If such states have also larger constant part of the diabatic energy 
a m > «o they can grow in the first half of the contour as exp(C'(4>)a m \t\), i.e faster then the 
initially filled state ipo, but being initially vanishing, amplitudes of such states can grow only 
due to transitions from the other states. At first half of the time-contour they are coupled 
only to states that are suppressed by much stronger exponents exp(C(0)/3j|t| 2 /2), < 0); 
therefore we do not expect that they become large in comparison with ipo up to = tt/2. 
In the second part of the path tt/2 < < states with such an asymptotic exp(— iat) 
already decrease exponentially and become suppressed in comparison with ipo; therefore we 
can expect that they do not break the approximation (jlj) for the state ipQ and have vanishing 
amplitudes at the end of the evolution. This is exactly in agreement with (J3J). 

Our arguments in support of (J2J) and (J3J) are certainly very intuitive and every step in 
the mathematically rigorous proof requires more detailed justification. However, we note 
that is also confirmed by all known exactly solvable classes featuring the possibility of 
counterintuitive transitions, namely by the Demkov-Osherov model^, the generalized bow- 
tie model^ and the model of two crossing bands of parallel levels^. Besides, we performed 
a number of numerical simulations with arbitrary choices of parameters. As we found, all 
they support our hypothesis ©• For example, in Fig|3]we show the time-dependence of the 
probabilities to find the system at states 2 and 3 in the model demonstrated scematically 
in FigE] if initially only the state 1 is occupied. One can deduce that generally during the 
evolution these probabilities can be rather high (> 0.1) and show oscillating behavior, but 
asymptotically at t — > +oo they vanish. Numerically we can simulate the evolution only in 
the finite time interval. For the evolution from t = —500 to t = 500 and the same parameters 
as in FigHJwe find |S 2 i| 2 = 5.18xl0- 7 and |S 31 | 2 = 3.11xHr 7 . In comparison \S U \ 2 = 0.234, 
15*41 1 2 = 0.295 and | -S' 51 1 2 = 0.472. We also note that although counterintuitive transitions 
have vanishing probabilities, the presence of the states 2 and 3 does affect other elements 
of the scattering matrix. Thus if we set all couplings of states 2 and 3 with all other states 
to zero, then numerically calculated nondiagonal transition probabilities are |S4i| 2 = 0.672 
and \Sz,i\ 2 = 0.094 that is different from our previous numerical result. 

As another example, consider a 4-state model shown in Fig0J Obviously, for e > the 
transition from the state 1 to the state 2 is counterintuitive but for e < it is not. Fig IS] 
shows numerically calculated final probabilities to find the system in all 4 states for the 
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FIG. 3: Time dependence of the counterintuitive transition probabilities for the model in Fig^ 
Triangles correspond to P(t) = |S , 2i(t, — oo)| 2 and boxes show P(t) = \S 3 i(t, — oo)| 2 . The choice 
of nondiagonal elements of the Hamilton operator is H12 = His = H23 = 0, H34 = 0.8, ^35 = 
0.3 + 0.24*, H 2A = 0.1 + 0.7*, H 25 = 0.5 + 0.1*, H 14 = 0.4+0.12*, H 15 = 0.25 + 0.2*, H Ab = 0.6 + 0.9i. 
The other elements are obtained by employing Hermitian properties of the matrix H. 
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FIG. 4: Diabatic levels of a 4-state model. The matrix elements of the Hamilton operator are 
chosen as follows. H\\ = -t, H 2 2 = —t - e, H 33 = t, #44 = 0M - 0.5, Hi 2 = 0, H 13 = 0.4 - 0.1*', 
H14 = 0.6, H23 = 0.4 + 0.5i, H24, = 0.2 + 0.3*. The other elements are obtained by employing 
Hermitian properties of the matrix H. 

evolution from t = —600 to t — 600 when initially only the state 1 is populated. 

One can see that the probability to remain in the state 1 does not depend on e, in 
agreement with the Brundobler-Elser conjecture. A tiny diviation from the Brundobler- 
Elser formula can be seen for two points with e closest to zero. However, this should be 
explained as due to the fact that e = is the critical point and it takes much more time 
for probabilities to saturate in its vicinity, but in simulations the time interval had to be 
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FIG. 5: Transition probabilities to all states in the model in Fig0]as functions of the distance e 
between levels 2 and 1. Long before level crossings the probability of the state 1 is set to unity. 
The transition probabilities are represented correspondingly by boxes for the state 1, stars for the 
state 2, crossed lines for the state 3 and triangles for the state 4. 

finite. At negative e all other probabilities strongly depend on the distance between states 
1 and 2. This can be explained partly even in the independent crossing approximation 
as due to the interference among different semiclassical paths leading to the same final 
state. However, for e > the independent crossing approximation does not predict any 
dependence of probabilities on e if initially only the state 1 is populated. Nevertheless, one 
can deduce from the FigElthat in addition to the state 1 only the transition probabilities to 
the state 2 become flat and have indistinguishable from zero magnitudes in agreement with 
(J3J). Transition probabilities to states 3 and 4 strongly depend on e there. This indicates 
that simulations were performed for the range of parameters where the independent crossing 
approximation fails unless its predictions become for some reasons exact. 

In conclusion, the generalization of the Brundobler-Elser hypothesis is proposed that 
states that counterintuitive transitions in the multistate Landau-Zener model with the linear 
time-dependence of diabatic energies are asymptotically forbidden. It is confirmed by all 
numerical tests and by all known exact solutions. In addition, we demonstrated that this 
result can be explained by continuation of the time-path into the complex plain, i.e. by 
the same approach as the one proposed by Landau to solve the two state model. As in any 
known exact solution of the multistate Landau-Zener model, the formula (j2J) coincides with 
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predictions of the independent crossing approximation. This fact points to the common 
origin of all such exact results. 
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